Energy Anomaly and Polarizability of Carbon Nanotubes 
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The energy of electron Fermi sea perturbed by external potential, represented as energy anomaly 
which accounts for the contribution of the deep-lying states, is analyzed for massive d — 1 + 1 Dirac 
fermions on a circle. The anomaly is a universal function of the applied field, and is related to 
known field-theoretic anomalies. We express transverse polarizability of Carbon nanotubes via the 
anomaly, in a way which exhibits the universality and scale-invariance of the response dominated by 
7r electrons and qualitatively different from that of dielectric and conducting shells. Electron band 
transformation in a strong-field effect regime is predicted. 

PACS numbers: 71.10.Pm, 85.35.Kt, 64.70.Rh 
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Fermion anomalies are universal contributions to low- 
energy properties of field theory originating from the bot- 
tom of the filled Fermi-Dirac sea. The primary examples 
in high-energy physics are the spontaneously generated 
photon mass in d = 1 + 1 , and the Adler-Bell- Jackiw 
chiral anomaly relating the decay of pion into two pho- 
tons to the number of quark colors Q . Several instances 
of anomalies are known in solid statephysics, with the 
chiral anomaly manifest in transport |4j and in fermion 
number fractionalization jjj, and the parity anomaly in 
d = 2 + 1 linked to the quantum Hall effect Q . 

Here we describe a new manifestation of fermion 
anomaly, appearing in the total energy of Fermi sea in 
the presence of an external field. While the net energy 
of Fermi-Dirac vacuum depends on ultraviolet cutoff, i.e. 
on the behavior at the bottom of the band, we find the 
external field-dependent part of the energy to be cutoff- 
insensitive. The latter energy is naturally divided into 
two parts, one given by a sum of the energy shifts of 
filled states near the Fermi level, and another, equally 
important, desribing cumulative effect of the states deep 
below the Fermi level. The latter, anomalous part, traced 
to Schwinger anomaly 0, has a universal form and can 
be expressed through the properties near the Fermi level. 

The Fermi sea energy field-dependence is the physical 
quantity central for many physical properties of materi- 
als. As an application, below we consider the response 
of the carbon nanotube (NT) n electron band to perpen- 
dicular electric field. The aforementioned separation into 
the normal and anomalous parts enables one to handle 
energy in a fully general and, at the same time, case- 
specific way, taking full account of the level quantiza- 
tion, curved geometry, spatial inhomogeniety, etc. The 
NT 7r-band is described by a tight-binding model on a 
honeycomb lattice Q. Near the Fermi level, at the tt- 
band center, the NT electrons are described by 2 + f 
Dirac fermions moving on a cylindrical surface |8j ■ Dirac 
model provides a simple analytic description of NT cur- 
vature and chirality 0, 0] , and of the effects of external 
fields [lO] in both semiconducting and metallic NTs. 



It might seem that understanding properties such as 
NT polarizability should require a detailed knowledge 
of the 7r electrons behavior on the lattice constant scale 
[nj . Here we find that, to the contrary, the energy cal- 
culated from Dirac model, with the anomaly properly 
accounted for, not only is numerically accurate, but also 
provides new insight into NT properties. We identify 
the dominant role of it electrons response compared to 
that of other orbitals, and explain the origin of the scale- 
invariant depolarization, independent of the NT diameter 
[ill H2I fl3j | , qualitatively different from that of metallic 
and dielectric shells. Our approach, with electron inter- 
actions included in an RPA fashion [lj] , is not limited to 
linear response: we apply it to study NT electron band 
transformation in the field-effect regime. 

The origin of the energy anomaly is exhibited most 
clearly by the example of chiral Id fermions on a circle 
[0,27r], described by the Hamiltonian H — —ide + U, 
U = 2a cos 0. In the Fourier representation \n) — e m6 , H 
is given by an infinite three-diagonal matrix: H nn = n, 
7i nn ±i — a. Eliminating the potential U by a gauge 



transformation ipiff) = e 



-2ia sin 



ijj(0), H = -id e , 



that the eigenvalues of TL are integers independent of U. 

The energy anomaly arises when the interaction U is 
truncated at a certain energy scale. To that end, let us 
consider a more general three-diagonal matrix 

W„ n = n, TL nn+ i —Hn+in = a n , n e Z, (1) 

with the sequence a n having different limits at n — > ±00: 
&n— »-oo = 5 a n~++oo — Although now the energy 
levels depend on a n , the above argument for spectrum 
robustness at constant a n indicates that this dependence 
is exponentially small at large |n|. The level shifts are sig- 
nificant for only a relatively small cluster of levels around 
n ~ n» where the switching of a n from to a occurs. 

Notably, the sum of all level shifts, StrH = ^ n Se n , 
depends only on the asymptotic values a n ->± 00 , while 
other details of the sequence a n do not matter. To see 
this, we truncate the matrix H at some large positive 
and negative n = ±iV, in which case the trace of this 
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(2N + 1) x (27V + 1) matrix is finite and explicitly a- 
independent. Since at JV > n„ the mutual influence of 
the levels at n ~ n* and n ~ N is exponentially small, the 
effect on tr?i due to truncating at n ~ TV is negative of 
that due to a n switching at n ~ n* , both being universal. 
[The levels at n ~ —A are unaffected by the truncation 
since 0n->— oa — and TL is diagonal.] Thus the sum of 
the level shifts at n ~ n* depends only on a n= +oo = a, 
giving a cutoff-independent contribution to StxH. 

The universal value E anom = StrH can be evaluated 
using slowly varying a n , \da n /dn\ <C |o«|. In this case, 
since the levels are unperturbed by constant a n , the level 
shifts will be small, which warrants using pertubation 
theory. Gradient expansion of a n in the vicinity of n — 
yields H = -id e + 2a cos 6 + b{e ie (~id e ) +h.c.}, with 
a = a n „ — bn*, b = da n /dn\ n=nr . The gauge transforma- 



tion ip(6) = e 



-Ha sin 9 



ip(0) transforms H into 



n 



id g + b {e i9 (-id e - 2acos6) + h.c.} . (2) 



The energies e n with n near n* obtained in the lowest 
order of perturbation theory are e n = {n\H\n) = n — 2ba. 
The sum of these level shifts, given by a full derivative, 
depends only on the asymptotic of a n : 



dn 2a, 



dOn 

dn 



(3) 



The relation of this result with fermionic energy emerges 
when one considers quenching of the external field a n = a 
coupling to the states at the Fermi sea bottom, modeled 
by a n < Ut — 0. We find that the anomalous contribution 
-E-anom is universal, i.e. it depends only on the properties 
near Fermi level |15| . The anomaly contributes additively 
to the energy along with the contributions due to fermion 
mass and confinement (see below). 

Electrons in a nanotube, a cylinder of radius R, are 
described by 2 + 1 massless Dirac model || . The states 
in a transverse field, labelled by momentum k along the 
tube, can be viewed as massive 1 + 1 Dirac fermions on 
NT circumference, the circle < y < 2irR: 



-ifrvady + Tivkfi + U{y) . 



(4) 



with U (y) the transverse field potential and Tivk the Dirac 
"mass". (Here a = —cr x and = a y 0.) We assume 
generic quasiperiodic boundary conditions ip(y + 2nR) — 
e 27TtS ip(y). At U = 0, the energy levels 



e£ = ±^J A R (n + S) 2 + (Hvk) 2 , Ar = hv/R, 

describe NT minibands. The phase 5 determines NT 
properties: S = ±| for semiconducting NT, (5 = for 
metallic NT, and \6\ <C 1 for the tubes with small gap 
induced by curvature 0] or by parallel magnetic field |lOj |. 

For the Fermi sea energy at finite U, the naive an- 
swer would be the sum over occupied states W — Y^, ^ e 7i i 



5e n = e n |j/ — e„|;7 = o- The level shifts Se n become 
very small away from the band center, at |n| 3> kR: 
In this limit, with the contribution of finite mass be- 
ing negligible, the Dirac problem (@J decouples into two 
chiral fermion modes, each having [/-independent spec- 
trum. However, despite the absence of level shifts, the 
states with large n contribute to the [/-dependent energy 
via anomaly due to the bandwidth cutoff, (hv/R)n* ~ 
— 10 eV for Carbon. As described above, the anomaly 
depends only on the properties near the band center: 



N 



f 



2ttHv 



2nR 



U 2 (y)dy, 



(5) 



where Nf = 4 is the number of electron flavors asso- 
ciated with spin and the points K, K' 7]. Note that 
the energy is additive for multiple fermion flavors, 
contrary to fermion-doubling cancellation typical of the 
chiral anomaly effects [lij . 

One can interpret Eq. JSJ , somewhat loosely, as a coun- 
terterm which eliminates the nonphysical contribution of 
infinite Fermi sea in the model i.e. the states out- 
side the Carbon band. A photon mass m 2 = e 2 /ir ap- 
pears in d = 1 + 1 QED l] after integration over massless 
fermions. Eq.JSJ exhibits a similar effect in the Dirac sys- 
tem in an external electromagnetic field eA^ = (U(jj), 0), 
with the mass term Jd 2 x ^m 2 ^ = — J E anom dt in the 
action. [We point out a distinction of our approach and 
Coleman's analysis 17] of massive Schwinger model in 
external field, which yields an effective action containing 
the field intensity rather than potential.] 

The energy anomaly is related to the chiral anomaly 
in d = 1 + 1, since in this case the chirality 7 5 = a enters 
the Hamiltonian Formally this can be seen via the 
change in the functional measure [l^. First note that a 
chiral gauge transformation [h = v = 1] 



tl>(y) = e-*'+W${y) , </> = / dy'U(y') 



(6) 



preserves the boundary conditions on tp, turning (J1J into 



'Hd 



-iady + ke 2i&4,{v) (3 . 



(7) 



Consider now the variation of the background field U (y) 
by 5U(y), related to the infinitesimal transformation 
p ; Tjj> = e *Tt 5 Hv)ip, \( y ) = J v dy'5U(y'). The corre- 
sponding Jacobian J = exp{2i J dtdy AX(y)} ^3 changes 
the action by —i In J = -8E anom t, yielding SE anom = 
i JdyXdyU, with A = tr 7 5 = = ~d y U(y)/2n 

the 2d anomaly Integration by parts gives the vari- 
ation 5E a „ nrrj = — i Jdy U5U, which is equivalent to the 
result © |23 . This derivation can be extended to include 
electron interactions, proving robustness of the anomaly. 
Turning to the NT response to transverse electric field, 



U(y) = -e£Rcos(y/R), 



(8) 
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we relate polarizability to the sum of Stark shifts 



From bottom to top: u = 0.2, 0.4, ... , 1.8 



E (k)= J2 5e ^ k ) 



(9) 



n— — oo 



6e n (k) = e n (k)\s — £n(fc)|£=o, taken over all occupied 
states, with k the electron momentum along the tube. 
We obtain the shifts 8e n (k) from the Hamiltonian 
In this calculation, the anomaly J5J) must be added to 
account for the finite band cutoff, formally absent in Q . 

The main effect of electron interaction is depolariza- 
tion, i.e. screening of the field inside NT. To obtain the 
RPA screening function 0, 0] of NT cylinder, we first 
show how the problem is reduced to the calculation of 
electron energy in the presence of an external field. The 
Gauss's law relates the fields inside and outside the tube 
with the induced surface charge density <r, 



£ cxt = £ + I • 4tt • N f a, <7 = V/(ttR 2 



(10) 



where V is the dipole moment per flavor and per unit 
NT length, and the factor 1/2 accounts for depolariza- 
tion in cylindrical geometry. In Ea. Hl(J|) we projected the 
actual charge density on the cost/? harmonic, (p = y/R, 
as cr(ip) — > Nf a cosip, ignoring higher order harmonics. 
The problem is then reduced to evaluating the dipole 
moment, given by V — —dW(£)/d£, where W(£) is one 
fermion flavor energy as a function of the inner field, 



W = NJ 1 



{E (k) + E a 



dk 
l} 2^ 



(11) 



Combining this with Ea. HlOfl . and using dimensionless 
u = e£R/An, we obtain u cxt = u + 2Nfj-V(u). Once 
the dipole moment V(u) is known, this relation deter- 
mines the inner field u in terms of the outer field u cxt . 

Taking the Stark shifts Se n (k) to the lowest order in £ 
and evaluating the integral over k (see Fig^l, we obtain 



W = 



0.196.. 
0.179.. 



for 6=1/3, 
for 6 = 0. 



(12) 



Notably, the linear relation V = au holds up to very 
large fields u ~ 1 (Fig^ inset), giving the depolarization 



(l + 2JV>«£)£. 



(13) 



With e 2 /hv = 2.7 this gives £ cxt /£ = 5.24 for 5 = 1/3, 
and £ cx t/£ — 4.87 for i5 = 0, in excellent agreement with 
the full tight-binding calculations |Tlj | . 

The screening (|13|) is radius-independent and is nearly 
the same in the metallic and semiconducting NTs. The 
latter is not surprising: screening is absent in a single 
Id mode approximation, sicne polarizability is related 
to dipolar transitions between different subbands. The 
scale-invariance of (|13fl obtained for a hollow NT cylinder 
resembles depolarization in a massive dielectric cylinder 
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FIG. 1: Partial dipole moment P(k) = —d(Eo(k)+E a , I1 om)/du 
as a function of k, obtained from lilt for semiconducting 
NT, where u = e£R 2 /hv is dimensionless field. Note the 
cancellation between Eo(k) and the anomaly ^ at kR ^> 
1, enforcing convergence of V = J P(k)dk/2n. Note also 
that P(k — > 0) is dominated by the anomaly, since Eo = 
at k = due to the chiral gauge invariance (|SJ. Inset: 
Dipole moment V per one fermion flavor versus u for metallic 
and semiconducting NT. Straight lines represent weak field 
linearization 11121 . Arrow marks the field u ~ 1.2 for which 
velocity changes sign in metallic NT (see text and Figl^J. 



|21|, |22j, |23j ■ The dipole moment of 7r electrons, found to 
scale with R 2 , should be contrasted to V cx R for hollow 
dielectric shell. The universal scale-invariant result l|13f) 
reflects the dominant role of 7r electrons in transverse 
response as compared to other Carbon orbitals. 

To emphasize the role of anomaly in this calculation 
we note that omitting this contribution would have led 
to a wrong sign of the response and also to a divergence. 
Indeed, due to an upward shift of the filled levels e~ 
(Fig|2J), we have Vq = dE /du > 0, corresponding to 
unphysical "diamagnetic" polarization sign. Also, the 
fc-dependence Eo(k) causes an ultraviolet divergence in 
the integral V = J P(k)dk/2-K, since E${k) increases with 
\k\, saturating at \k\R 3> 1 at an asymptotic value \u 2 . 
Both difficulties are resolved by taking into account the 



negative E a 



1„,2 



The resulting integral i|ll|l con- 



verges after E (k) is offset by -Eanom (Fig0. 

To illustrate the effect of transverse field, here we 
examine the NT electron spectrum. The changes are 
most dramatic in a strong field 0] which mixes differ- 
ent NT subbands, u ~ 1, or e£R ~ A R , £ [MV/cm] ~ 
5.26/R 2 [nm 2 ]. In metallic NTs the electron velocity 
v = de/dp decreases and can even reverse sign, causing 
Fermi surface breakup. This could lead to interesting 
many-body effects such as the Luttinger correlations in- 
crease due to enhanced e 2 /hv, or instability with respect 
to exciton formation for the negative-u states. Semicon- 
ducting NTs exhibit the effective mass sign change, ac- 
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FIG. 2: Electron bands transformation in the strong field- 
effect regime, u = e£R/AR ~ 1. Top; Velocity reversal in 
metallic NT occurs at u > u e ~ 1.2. Bottom: Effective mass 
sign change in semiconducting NT at u > u c w 0.62. The 
bands are shown for u below and above the critical value. 
Inset: Energy gap suppression in a semiconducting NT. 

companied by field-induced suppression of the excitation 
gap (Fig|2J) that could be manifest in activated transport. 
The above screening calculation provides an estimate for 
the fields necessary to observe these effects. The rel- 
atively small outer-to-inner field ratio ~ 5, as well as 
availability of NTs of a few nm radius, puts the required 
voltage across the tube in a feasible range of a few Volts. 

In summary, we considered the energy of Fermi sea 
for Dirac fermions in an external field, expressed via en- 
ergy anomaly through low-energy properties. The en- 
ergy anomaly, applied to the polarizability of nanotubes, 
provides insight in their response properties, notably 
the scale-invariance of depolarization, the difference of 
screening in semiconducting and metallic NT, and rela- 
tive importance of the 7r band compared to other Carbon 
bands. Electron bands exhibit dramatic change in the 
strong field-effect regime. 
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